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Anderson models, based on quantum chemical studies of 
the molecule of U(CsH8)2, are applied to investigate the prob- 
lem of an U impurity in a metal. The special point here is 
that the U 5/-orbitals are divided into two subsets: an al- 
most completely localized set and a considerably delocalized 
one. Due to the crystal field, both localized and delocalized U 
5/-orbitals affect the low-energy physics. A numerical renor- 
malization group study shows that every fixed point is char- 
acterized by a residual local spin and a phase shift. The latter 
changes between and n/2, which indicates the competition 
between two different fixed points. Such a competition be- 
tween the different local spins at the fixed points reflects itself 
in the impurity magnetic susceptibility at high temperatures. 
These different features cannot be obtained if the special char- 
acters of U 5/-orbitals are neglected. 



In the field of heavy- fermion systems, the research area 
has been extended from rare-earth systems (mainly Ce 
based compounds) to actinide ones (mainly U). From the 
numerous studies of Ce systems Q , we have learned that 
the on-site physics (Kondo problem) is essential for un- 
derstanding the formation of heavy quasiparticle bands. 
In the case of a Ce 3+ ion, for which the 4/ 1 configu- 
ration is most stable, a singlet state may be formed at 
the Ce site via the hybridization between the conduction 
band and the f-orbitals. This helps us to understand one 
possible origin of heavy quasiparticles in lattice systems, 
although other mechanisms have been identified in dif- 
ferent materials such as YD4AS3 and Nd2- a; Ce 2 ,Cu04 

On the other hand, in the case of an U, we cannot 
easily extend the on-site physics to the lattice case be- 
cause of the complicated atomic structure. In real ma- 
terials, U has a 5/ 2 or 5/ 3 configuration. In order to 
understand the atomic physics of the U site, we must 
take into account the Hund's coupling as well as direct 
Coulomb interaction. As for Ce, the effect of the crys- 
tal field is important for U; however, it plays a differ- 
ent role here. The crystal field deforms the degenerate 
5 /-orbitals and brings anisotropy to the hybridizations 
between the deformed orbitals and conduction electrons. 
Kusunose and Miyake studied how the Hund's coupling 
affects the Kondo effect in such impurity systems [Q. 
However, they confined themselves to two-fold degen- 
erate orbitals in the U 5/ 2 configuration: one of them 
shows a hybridization, while the other is completely lo- 



calized. Their Kondo model, to which the numerical- 
renormalization-group (NRG) method §-§] was applied, 
leads to the result that the effect of Hund's coupling- 
is irrelevant and the fixed point is of the strong cou- 
pling type. The irrelevant coupling between a local spin 
and conduction electrons is ferromagnetic near the fixed 
point. Such anisotropic systems were also studied by 
others in the limit of strong spin-orbit coupling, where 
the two-channel Kondo effect was discussed Jj|§]. Al- 
though we have to learn the atomic structure of U in 
metals to present a more realistic model, it has not been 
specified yet on the experimental side. However, some 
information which is useful to solid state physics can 
be extracted from accurate quantum chemical calcula- 
tions of some molecules with relatively large ligands, e.g., 
Ce(C 8 H 8 ) 2 go| and U(C 8 H 8 ) 2 Q. The ground state of 
Ce(C 8 H 8 ) 2 was found to be an open-shell singlet i}A lg ) 
JTo| , which is analogous to a Kondo singlet of magnetic 
impurities in metals. A splitting between the lowest sin- 
glet and excited triplet states corresponds to a low-energy 
scale for heavy fermion systems in Ce cases. According 
to recent results based on the multiconfiguration self- 
consistent field (MCSCF) calculations of U(C 8 H 8 ) 2 [0, 
the U 5/-orbitals can be divided into two subsets. Two 
of the seven U 5f orbitals, which have the same symmetry 
(e 2u in T) S h point group) as the ligand HOMO and thus 
strongly hybridize with the ligands, are significantly de- 
localized (~ 90% /-character). The other five 5f-orbitals 
have different symmetries (ei„, e^ u , a 2u ) from that of the 
ligand HOMO, and are thus almost completely localized 
(> 99% /-character). This extent of delocalized char- 
acter is not found for lanthanocenes such as Nd(C 8 H 8 ) 2 , 
where all the 4/-orbitals are almost completely localized. 
First, we design a model Hamiltonian to reproduce the 
low-lying states of U(C 8 H 8 ) 2 in order to understand 
the important aspects. 

Table I (a) shows a few low-lying states of U(C 8 H 8 ) 2 
taken from ref. |ll[]. The low-lying states can be repro- 
duced by using an Anderson model (Q) where the different 
characters of the 5/-orbitals are taken into account: 



H = Hf + Hi, 
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Here the energy is measured from the ligand orbital. Hf 
represents the energy of U 5 /-electrons; Hi is the hy- 
bridization Hamiltonian between the 5/ and ligand or- 
bitals. fwMj (fmcr) and l\_ 2 a {l±2,a) are the creation (an- 
nihilation) operators for the 5/ and ligand electrons with 
spin a, respectively. The orbitals are denoted by the 
magnetic quantum number m (= 0, ±1, ±2, ±3 for /). 
Only one pair of orbitals with m — ±2 has a noticeable 
hybridization jOJ. The second and third terms in Hf 
are the direct and exchange interactions among the 5/- 
electrons, respectively. In this Hamiltonian, every term 
depends on the orbitals deformed by the crystal field in 
the molecule. Although the spin-orbit coupling, which is 
neglected here, is relatively large for heavy atoms like U, 
it was found to be of minor importance compared with 
the Hund's coupling. As far as we assess the low-lying 
energy levels, it only yields a slight modification [flTfl . 

The low-lying states shown in Table I (a) can be re- 
produced qualitatively with our Anderson Hamiltonian 
where the parameters are derived directly from quantum 
chemical calculations. Since there are too many param- 
eters in the Hamiltonian ([!]), we consider two kinds of 
parameter sets. In both sets, Eq, e±i, and £±3 are lower 
than e±2- In Set I, a fine anisotropy is introduced to the 
localized orbital energies, which results in a small split- 
ting of the low-lying states: 

£0 = -1.205, £±i = -1.200, 

£ ±2 = -1.150, £ ±3 = -1.210, 

U ,a = U Xtl = U 3 , 3 = C/3,-3 = U QA = 0.700, 

Uo,3 = U± t -i = Ui,3 = C/ 1,-3 = U() t 2 = U\.2 = U\-2 

= U 2 , 3 = U 2 ,- 3 = U2.2 = U2.-2 = 0.650, 
2J m , m > = 0.040, V± 2 = 0.032, (4) 

where the values are measured in atomic units (au= 27.2 
cV). The Coulomb couplings satisfy 



Urrt.rr).' — Um.'m, — U— 



(5) 



In Set II, the localized orbitals are degenerate, while more 
precise values derived from ab initio HF calculations are 
given to the Coulomb couplings: 

£0 = -1.200, £±i = -1.200, 
£ ±2 = -1.150, £ ±3 = -1.200, 

U ,o = E7i,i = C/3,3 = 0.700, C/3,-3 = 0.71, U ,x = 0.690, 
U ,3 = U h3 = U 1; - 3 = U x>2 = C/i,_ 2 = U 2 , 2 = 0.650, 
U 2 ,3 = U 2 ,- 3 = U 2 ,-2 = 0.640, C/i,_i = 0.660, 
U , 2 = 0.630, 2J m>m , = 0.040, V± 2 = 0.032. (6) 

If the fine anisotropy is not taken into account, the 
ground state is still 21-fold degenerate. This is not the 
case in Table I (a). We note that the first and second 
excited states derived from the Anderson model are de- 
generate (Table I (b)), which are actually split due to 
D sh symmetry in U(C 8 H 8 ) 2 . 



We then can present an Anderson model to study an 
U impurity in metals, replacing the ligand electrons by 
conduction electrons. In order to bring out the orbital 
effect, we reduce the number of orbitals: only one de- 
localized and one two-fold degenerate localized orbitals 
are retained. The delocalized orbital, which is denoted 
by M — 0, represents the orbitals with m = ±2. In 
the same way, the localized orbitals, which are denoted 
by M = ±1, represent those with m = 0, ±1,±3. 
Anisotropics in the Coulomb and Hund's couplings (see 
(|]) and (H)) do not affect the low-energy physics and do 
not change the high-energy physics qualitatively. The 
most important quantity is the splitting between the lo- 
calized and delocalized orbitals. Then the simplified An- 
derson Hamiltonian is written as 



H — Hk + Hf + #hyb, 
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where Hy represents the energy of conduction electrons. 
a \ a ( a k,o.o-) is the creation (annihilation) operator for 
the conduction electrons with wave number k, orbital 
quantum number M = and spin a. 

The Anderson Hamiltonian (fjj) is transformed to a 
hopping type of Hamiltonian via the standard procedure 
of the NRG theory The transformed Hamiltonian 
satisfies a recursion relation of the form 

H N+1 = A 1/2 H N + XX s ko,<r s ^+iA<x + h.c), (11) 

a 

and the impurity part is given by 



H = 
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where the energies are measured in units of the half width 
of the conduction band D. a {s„.o y a) corresponds to 
the creation (annihilation) operator for the conduction 
electrons. For the logarithmic discretization parameter, 
A = 2 ~ 3 is used. We will fix the strength of hybridiza- 
tion, the Coulomb and Hund's coupling at T = 0.100, 
U = 0.650 and Jh = 0.020, respectively, throughout the 
remaining part of this letter. Although these values are 
not connected directly with those given in (0) and (pf) , the 
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ratio of Coulomb and Hund's couplings is based on the 
analysis of the model (|l|) . For practical reasons, the hy- 
bridization is assumed to be relatively high, which should 
not change the low-energy physics as compared with a 
smaller hybridization. 

As a result of our NRG calculation, several types of 
fixed points are obtained, which depend sensitively on 
the energy splitting between localized and delocalized or- 
bitals. Every fixed point can be characterized by a resid- 
ual local spin (in Fig. 1) and a phase shift (in Fig. 2). 
The latter will be explained later. Figure 1 shows the 
size of the local spin at the fixed points depending on the 
choice of parameters. The f 2 configuration is found to 
be more stable than any other one within the parameter 
range -2(U - J H ) = -1.26 < e < -{U - J H ) = -0.63, 
provided £*o — £±i- It is most stable around a point 
e = e±i = —0.945. When both orbital energies are far 
from this point, charge fluctuations increase between f 2 
and f 1 or f 2 and / 3 . When eo > —0.63, f 1 is more stable 
than f 2 , while in the region e~o < —1.26, we find that f 3 
is more stable than / . 

We can reproduce the excitations near the fixed point 
by using the following effective Hamiltonian: 



H 
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with a ferromagnetic exchange interaction between the 
residual local spin and the quasiparticles defined at the 
fixed point. The positive exchange coupling w(N) de- 
creases as N increases. This type of Hamiltonian was 
studied by Kusunose and Miyake in the case of the under- 
screened Kondo effect j| . Their discussion was restricted 
not only to the Kondo regime where charge fluctuations 
are eliminated but also to the special case of our model 
where e"o = £±i is satisfied. In general, the second term in 
( |l3| ) is marginal, while w(N) vanishes at the fixed point. 
The excitations from the ground state can be described 
by a single-particle Hamiltonian 



H 



N 



Vi h]hi] 



(14) 



where rw ) is the /-th excitation energy of a particle 
(hole). The latter is given by 



A i-i T 5A (jv : odd) 
A}- 1 / 2 **/* (A: even) 



(15) 



where 1 -C I *C N. The phase shift 5 depends on e*, 
namely, e~o and e±%. As shown in Fig. 2, 5 changes 
monotonously as far as the residual local spin keeps its 
size. This behavior depends only on e~o when T, U and Jh 
are fixed. In other words, the localized orbitals (M = ±1) 



are only relevant to the size of the local spin at the fixed 
point. Usually, <5 = implies that the local spin de- 
couples from the conduction electrons, while S = n/2 
means that the size of the local spin shrinks from S to 
S — 1/2 due to a spin compensation. The change of 
phase shift between and it/2 indicates a competition 
between the two fixed points [ p^[ . The maximum and 
minimum of S depend on A — e±\ — eo, which changes 
the boundaries of the phase diagram in Fig. 1. In the 
vicinity of the boundary in Fig. 1, such competition is 
observed clearly in the impurity magnetic susceptibility 
Ximp at high temperatures. In Fig. 3, a shoulder appears 
in Txi,„p within the range 0.001 < T/D < 0.01 (T is 
temperature). At lower temperatures, the susceptibility 
shows a Curie law: Txi mp = 2/3 for S = 1 at the fixed 
point, while Tximp — 1/4 when S = 1/2. 

Throughout this letter, much attention is paid to the 
important individual characters of U 5/ 2 orbitals: two 
of the seven orbitals are delocalized and the others are 
almost completely localized. Several low-lying states of 
U(CsHg)2 can be reproduced qualitatively by using a re- 
alistic Anderson model including such U 5/-characters. 
Our NRG study of the impurity problem based on a sim- 
plified Anderson model shows that every fixed point is 
characterized by a residual local spin and a phase shift, 
which depend sensitively on the splitting between the lo- 
calized and delocalized orbitals. The change of phase 
shift between 5 — and 6 — ir/2 is due to a competi- 
tion between two fixed points corresponding to the two 
limiting cases. For the 5/ 2 configuration, as a result of 
competition between 5 = 1/2 and S = 1 at the fixed 
points, a shoulder can be observed in Tx; mp at high tem- 
peratures. Our results in this letter could be extended 
to infinite dimensional lattices |l3| as was done for the 
Anderson model in the absence of orbital degeneracy [0 . 
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FIG. 1. Phase diagram of the residual local spins S at the 
fixed points: A = s±i — io (io < 0) is the splitting energy 
between the localized and delocalized orbitals. 

FIG. 2. Dependence of the phase shift 8 on the delocalized 
orbital energy io. Here i±i — io < is assumed. S changes 
abruptly when eo crosses the boundary of the phase diagram 
shown in Fig. 1: in this case, eo = —1.02. 

FIG. 3. Temperature dependence of the impurity mag- 
netic susceptibility. The localized orbital energy is fixed at 
e ±1 = -0.945. 
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